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1. Introduction 


In such excellent treatises on hydrodynamics as those of LAMB, [6], and BASSET, 
[2], we find that even when friction is neglected, long waves in a rectangular canal 
must necessarily change their form as they advance, becoming steeper in front 
and less steep behind'. Yet since the investigations of de BOUSSINESQ ({3]), Lord 
RAYLEIGH ((8]), and St. VENANT ([9]) on the solitary wave, there has been 
some cause to doubt the truth of this assertion. Indeed, if the reasons adduced 
were really decisive, it is difficult to see why the solitary wave should make an 
exception”; but even Lord RAYLEIGH and MCCowaAN ([7|), who have successfully 
and thoroughly treated the theory of this wave, do not directly contradict the 
statement in question. They are, as it seems to us, inclined to the opinion that the 
solitary wave is only stationary to a certain approximation. It is the desire to settle 
this question definitively which has led us into the somewhat tedious calculations 
which are to be found at the end of our paper. We believe, indeed, that from 
them the conclusion may be drawn, that in a frictionless liquid there may exist 
absolutely stationary waves and that the form of their surface and the motion of 
the liquid below it may be expressed by means of rapidly convergent series. But, in 
order that these lengthy calculations might not obscure other results, which were 
obtained in a less elaborate way, we have postponed them to the last part of our 
paper. First, then, we investigate the deformation of a system of waves of arbitrary 
shape but moving in one direction only, i.e. we consider one of the two systems 
of waves, starting in opposite directions in consequence of any disturbance, after 
their complete separation from each other. By adding to the motion of the fluid a 
uniform motion with velocity equal and opposite to the velocity of propagation of 
the waves, we may reduce the surface of such a system to approximate, but not 
perfect, rest. 

If, then, 1+ 7 (7 being a small quantity) represent the elevation of the surface 
above the bottom at a horizontal distance x from the origin of coordinates, we 
have succeeded in deducing the equation 
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Tt seems that this opinion was expressed for the first time by Arry, [1] 

Though the theory of the solitary wave is duly discussed in the treatise of BASSET, the 
inconsistency of his result with the doctrine of the necessary change of form of long waves seems 
not to have sufficiently attracted the attention of the author. 
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where qa is a small but arbitrary constant, which is in close connexion with the 
exact velocity of the uniform motion given to the liquid, and where 0 = ee a 
depends upon the depth / of the liquid, upon the capillary tension T' at its cuitane: 
and upon its density p. 

On assuming ou = 0 we of course obtain the differential equation for stationary 
waves, and it is easily shown that the well-known equation 


n =h-sech? x ui 
4o 

of the solitary wave is included as a particular case in the general solution of this 
equation®. But, in referring to this kind of wave, we have to notice the result that, 
taking capillarity into account, a negative wave will become the stationary one, 
when the depth of the liquid is small enough. On proceeding then to the general 
solution, a new type of long stationary wave is detected, the shape of the surface 
being determined by the equation 


h+k h 
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where fh is a constant and k with 0 < k < 1 is the elliptic modulus of the function 
cn, known as cosinus amplitudinis. 

We propose to attach to this type of wave the name of cnoidal waves (in analogy 
with sinusoidal waves). For k = 0‘ they become identical with the solitary wave. 
For large values of k they bear more and more resemblance to sinusoidal waves, 
though their general aspect differs in this respect, that their elevations are narrower 
than their hollows; at least when the liquid is not too shallow, in which latter case 
this peculiar feature is reversed by the influence of capillarity. For very large values 
of k these cnoidal waves coincide with the train of oscillatory waves of unchanging 
shape discovered by STOKES (({13]), which therefore in the theory of long waves? 
constitutes a particular case of the cnoidal form. Indeed the equation ({12], page 


3The functions secans hyperbolicus or hyperbolic secant and cosecans hyperbolicus or hyperbolic 
cosecant are hyperbolic functions and are defined according to: 
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ev’ —e-* — sinha 
They are per definition the inverse of hyperbolic sine and hyperbolic cosine. 
4see Appendix B, Table B.1 


°>STOKES solution is more general in so far as it applies also to those eases wherein the depth of 
the liquid is moderate or large in respect to the wave-length. 
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210) obtained by STOKES, when written in our notation, becomes 


=e Qra 3h?2d? Arar 
OND Glee KO 
but, as Sir G. STOKES remarks, in order that the method of approximation adopted 


by him may be legitimate, Ah must be a small fraction. Now, when capillarity is 
neglected, the wave-length A of our cnoidal waves is equal to 


4K . VB 
/3-(h+k)| 


with K = K(k) being the complete elliptic integral of the first kind and therefore 


dh 16K2h 16 
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This is a small fraction only when M, the modulus, is small, but the cnoidal waves 
then resemble sinusoidal waves; and it is obvious that in this case the equation of 
their surface may be developed in a rapidly convergent FOURIER-series, of which 
Sir G. STOKES has given the first two terms. 

After some more discussion about these cnoidal waves, concerning their velocity of 
propagation and the motion of the particles of fluid below their surface, we proceed 
to a closer examination of the deformation of long waves. To this effect we apply 
the equation for om to various types of non-stationary waves, and it will appear 
that, though sinusoidal waves become steeper in front when advancing, other types 
of waves may behave otherwise. 
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In our investigations (in accordance with the method used by Lord RAYLEIGH, ([8], 
page 257, whose paper has been of great influence on our researches) we start from 
the supposition that the horizontal and vertical velocities u and v of the fluid may 
be expressed by rapidly convergent series of the form 


w=ftyfity ft 
v=ybity dot: 


where y represents the height of a particle above the bottom of the canal, and 
where f, fi,---,; $1, @2,.-. are functions of the horizontal distance x and time t. Of 
course the validity of this assumption must be proved later on by the fact that 
series of this description can be found satisfying all the conditions of the problem. 
From one of these conditions, viz., the incompressibility of the liquid, which is 
expressed by ou + oe = 0), we may deduce 
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and, moreover, if ¢ be the velocity potential and w the streamfunction: 
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which set of equations satisfies for the interior of the fluid all the conditions of the 
problem, whilst at the same time it is easy to see that for long waves these series are 
rapidly convergent. Indeed, for such waves the state of motion changes slowly with 
x, and therefore the successive differential quotients with respect to this variable of 
all functions referring, as f does, to the state of motion, must rapidly decrease. 
Passing now to the conditions at the boundary, let p; (a constant) be the atmo- 
spheric pressure, p the pressure at a point below the surface where the capillary 
forces cease to act, and T the surface-tension. We then have, distinguishing here 
and elsewhere by the suffix (;) those quantities which refer to the surface, 
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but according to a well-known equation of hydrodynamics’, 


P' d¢, 1 
ra (t) - a —=-(ui+vt) —om, 
therefore 
Pi dd 1 2 2 T Oy 
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'The BERNOULLI equation, which reads in the case of unsteady potential flow: 


valid also for unsteady—or time dependent—flows. Here a denotes the partial derivative of the 


velocity potential ¢ with respect to time t, and v = |V¢| is the flow speed. The function y(t) 
depends only on time and not on position in the fluid. 
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By differentiation with respect to x, Equation 2.5 may be written as 
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Moreover, a second equation must hold good at the surface, viz. 
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In order to satisfy Equation 2.6 and Equation 2.7 by the method of successive 
approximations, we put y; =!+7, f = q+, where | and q are supposed to 
be constants, and 7 and ( small functions depending upon x and t. Dealing then, 
with the fact that for long waves, whose wave-length is great in comparison with 
the depth of the canal, every new differentiation with respect to x gives rise to 
continually smaller quantities, these equations become as a first approximation: 


OB | Op | AO <6 
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and are satisfied by taking 
dy _ db 
di die? 
= 2 (7 +a), 


and 
go = Val, (2.8) 


where a is an arbitrary constant which we will suppose to be small. 
It is obvious that this solution coincides with the one usually given for the case of 
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long waves of arbitrary shape made stationary by attributing to the fluid a velocity 
equal and opposite to that of the waves, on the assumption that the velocity in a 
vertical direction (v) may be neglected and that the horizontal velocity (wu) may be 
considered uniform across each section of the canal. 

But, if we wish to proceed to a second approximation, we have to put 


f=q-P-(nt+a+y), (2.9) 


where ¥y is small compared with 7 and a. On substituting this in Equation 2.6 and 
Equation 2.7 and on writing out the result, rejecting all terms? which are small 
compared with any one of the remaining terms, we find respectively: 
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In eliminating i from these equations, we have at last 
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This very important equation, to which we shall have frequently to revert in the 
course of this paper, indicates the deformation of a system of waves of arbitrary 
shape, but moving in one direction only. Before applying it, we may point out the 
close connexion between the constant a, which may still be chosen arbitrarily, and 
the uniform velocity given to the fluid. Indeed it is easy to see from Equation 2.1 
and Equation 2.9 how a variation da of the constant a corresponds to a change 
dq = *6a in this velocity, but, on taking the variation of Equation 2.12 with 
respect to a, we obtain 
5.97 _% 5.90 __ 5, On 
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which equation may be easily verified geometrically. 
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For stationary waves, 22 must be zero. Therefore we have from Equation 2.12 


> dt 
1 2 1 oF 
a (Sr'4 an+ ro: 1) <0. 
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and by multiplication with 6d7 and further integration 


Fa) 2 
co + 6e1n + 7° + 2an*? +o- (=) =, (3.2) 


If now the fluid be undisturbed at infinity, and if 1 be taken equal to the depth 
which it has there, then Equation 3.1 and Equation 3.2 must be satisfied by 7 = 0, 


a = 0, and ay = 0. Therefore, in this case c,; and cy are equal to zero, and 
Equation 3.2 leads to 
On / iP (n + 2a) 
— = +4/ -——___.. 3.3 
Ox o oe) 


Here, before we can proceed, we have to discriminate between o positive and a 
negative. In the first case 2a is necessarily negative because on must be real for 
small values of 7. If then, we put it equal to —h, we have 


dn yi 
Ee eae at, ao) Ppa 
dx ee ; 


from which, supposing x to be zero for 7 = h, we easily obtain the well-known 
equation of the positive solitary wave, viz.: 


/ h 
2 
7) h-sech* x (3 ) 
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In the second case 2a must be positive. In putting it equal to h, and in 
substituting —7 for 7, we have from Equation 3.3 


/ 
1 
OE a tyf—t a h—1, 
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n= —n! = —h-sech* x -4/-—. 
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This is the equation of a negative solitary wave, and we are able now to draw 
the conclusion that whenever o is negative; that is whenever the depth of the 


or, by integration, 


liquid is less than a the stationary wave is a negative one. For water at 
20°C this limiting depth is equal to 0.47cm (T = 72 dyn/cm , g =981cm?s71, 
p = 0.998 B.A.U.). 

Now, for a further discussion of Equation 3.2, we drop the assumption that the 
fluid is undisturbed at infinity. If then / be taken equal to the smallest depth of 
the liquid, we must have vt = 0 for 7 = 0, and therefore in virtue of Equation 3.2 
cz = 0. On supposing then o positive!, c, must be negative in order that oa may 
be real for small positive values of 7, but then the equation 


n’ + 2an + 6c; = 0 (3.5) 


has a positive root h and a negative —k, and we may get from Equation 3.2 


OH = fbn. (hn): (h-+ 1). (3.6) 


By substitution in this Equation 3.6 of 7 = hcos? x and by integration, we find 


n=h-cn’s- ae (moan fp4) (3.7) 


which is the equation of a train of periodic waves whose wave-length increases when 
k, decreases. 


/ 


lWhen a is negative, let then / be equal to the greatest depth. On substituting ¢o = —o’, n = —n 
we have again c, negative, 


and finally 


where h and —k are the roots of the quadratic equation 7’? — 2an + 6c, = 0 
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For k = 0 this length becomes infinite, and the equation may be shown to coincide 
with Equation 3.4. The following Figure 3.1 represents such a train of stationary 
waves for the case in which k = ah, M = 0.8. 


ee ef a oe 
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Figure 3.1.: Train of stationary waves 
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4. Stationary Periodic Waves 
(Cnoidal Waves) 


Proceeding now to a further investigation of the waves determined by Equation 3.7, 
we calculate from Equation 2.10 and Equation 2.11 the value of y. From these 
equations we get 


Dye. My Om Pls oP Vy 
3 0p) On?” 


or, by integration 


fs) ed EN OF 
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where the constant of integration is rejected because its retention would only have 
had the effect of augmenting in Equation 2.9 the value of the arbitrary constant a. 
On substituting, then, f from Equation 2.9 in Equation 2.1 and Equation 2.2, 
observing that in virtue of Equation 3.1 


O7n il 1 
Faz = Tag” Or + 4an + 6er) =—55 (30° —2-(h—k)-9— hk) , 


these equations are replaced by 


vv fb for 0B (aa) 
hh — Sn *|} — is {(h- k)-n+ shh 5} “yt. (4.1) 
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(4.2) 


When & = 0, they determine the motion of the fluid for a solitary wave. 


In the first place we now will endeavour to calculate the velocity of propagation. 
For the solitary wave this is simple enough. If we consider that the liquid at infinity 
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is brought to rest when a uniform motion with a horizontal velocity 


—q=-vVql- (1 + i) (4.3) 


is added to the motion expressed by Equation 4.1 and Equation 4.2, it is clear that 
this velocity, with reversed sign, must be taken for the velocity of propagation of 
the solitary wave. 

But for a train of oscillatory waves Sir G. STOKES has shown ([{12] ,page 202) that 
various definitions of this velocity may be given, leading at the higher order of 
approximation to different values. It seemed to us most rational to define it as 
the velocity of propagation of the wave-form when the horizontal momentum of 
the liquid has been reduced to zero by the addition of a uniform motion. This 
definition corresponds to the second one of Sir G. STOKES. According to it, we 


have to solve the equation 
» [+n 
| a | (u—q)dy =0, (4.4) 
0 0 


where g denotes the velocity of propagation, and where 


_ 2K Ve 
\= (4.5) 


is equal to the wave-length. 
If, then, 


v= fnar=s o_. {(h+k)- B(K) — kK} 


=A-{ (nek) A? eh (4.6) 


denote the volume of a single wave reckoned from above its lowest point, we get 
from Equation 4.4, retaining only such terms as are of the first order compared 
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with 7, h, and k: 


_ fede fy udy — fo {val — V/9n- 3: (k-h)- 9} tna 


iG dx fo dy Puen +) dx 
koh 
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= a-(1-S5" 5) 
a (1. (k +h) ee) 2). (4.7) 
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On subtracting this velocity from that expressed by Equation 4.1, we obtain 


= t.(9-7) (4.8) 


and it is obvious at once that in this manner we have annulled the velocity of the 
particles for which 
ce 
1 ec 

This last equation has a simple geometrical meaning. It designates those particles 
E (Figure 3.1) whose height above the bottom of the channel is equal to the 
height where the surface of the liquid would stand when the waves were flattened. 
Therefore, for a first approximation we may say that the various particles of the 
fluid change the direction of their horizontal motion at the very moment when one 
of these points EF is passing over them. 


We now proceed to the calculation of the path of a single particle of fluid. Let 
Xp and yo denote the coordinates of such a particle at the origin of time, and 
wv =2+€', y =yo+C its coordinates at the time t, u’ and v’ its horizontal and 
vertical velocity at that time, / +1’ its elevation above the bottom, then we have 


t t 
/ / g / V 
~ dt =,/2- eee 
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t t / 
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Here 1’ is equal to the value of 7 for x = x’ + qt; and therefore we have dx = 
(u’ — q) dt, or to a first approximation dt = = dx = 74 dx; but then 


' ee a h eae ee 
eae pe ae —— dx 
x0 


Or, according to a well-known formula! (see [4], ch. VI, §187, page 143) 


,  (h+k)-A 2K - (xo + Vg - t) 2K x9 
¢= GE a. ( )-2-( i if (4.9) 


At the same time we have 


; i totValt Ay h 3 (h +k) 
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Of course, as all fluid particles with the same y describe congruent paths, these 
formulee may be simplified by supposing xp = 0. 


1Z(u) =u(1- ze ) —M? fi'snudu 
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5. Deformation of Non-Stationary 
waves 


In order to study the deformation of non-stationary waves, we will now apply our 
Equation 2.12 to various types of waves. 


5.1. Solitary waves 
As a first example we choose a solitary wave whose surface is given by 
n = h-sech” px (5.1) 


According to Equation 2.12 the deformation of this wave is expressed by 


dn 3qoph 5 9 

— =— - (4op* —h)- |—sech 

a i ( op h) sech* px 
2 (a+ 2cp") 9 3 
=. ___—__| - sech - tanh ; 2 
a dep=h) sech* px - tanh” px (5.2) 


But before we are able to draw any conclusion from this expression, it is necessary 
to separate the two parts of sue of which the first is due to a true change of form of 
the wave-surface, whilst the second may be attributed to a small advancing motion 
of the wave, which is left after the addition of the uniform motion with velocity 
qo = Vgl. To this effect we have still at our disposal the quantity a, whose close 
connexion with the uniform motion, which we have added in order to make the 
wave nearly stationary, has been indicated above. 

One of the best ways to obtain the desired separation is certainly to make stationary 
the highest point of the wave, and this is effected by fulfilling the condition 


2+ (a + 2cp”) = 3- (4cp? — h) , 


or 


3 
= 467" = =h: 
a= dap’ — sh; 
for in that case Equation 5.2 is simplified to 


dn 3qoph— 
dt l 


(dcp? — h) « sech” px - tanh? pr (5.3) 
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Figure 5.1.: Change of form of a wave from Equation 5.3 


and then, for x = 0, 


ea = gan 0 
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In discussing this Equation 5.3, we see at once that a solitary wave according 
to Equation 5.1 is stationary when h = 4op?; and this is in accordance with the 
Equation 3.4 of the stationary solitary wave which we have obtained above. When 
h > 4op’, the change of form of the wave, calculated from Equation 5.3, is shown 
by the dotted line in Figure 5.1. Here the wave becomes steeper in front!, whilst 
for h < 4op? the figure would show the opposite change of form, when, contrary to 
the opinion expressed by Airy and others, the wave becomes less steep in front and 
steeper behind. 

If, now, we take account of the fact that, as may easily be inferred from Equation 5.1, 
the wave-surface becomes steeper in proportion as p is increased, we are then 
justified in saying that a solitary wave which is steeper than the stationary one, 
corresponding to the same height, becomes less steep in front and steeper behind, 
but that its behaviour is exactly opposite when it is less steep than the stationary 
one. 


5.2. Cnoidal Waves 


Applying Equation 2.12 to the cnoidal wave, 


n=h-cn’ pz, (5.4) 


'The left side of the figure is the front side of the wave, because the wave has been made stationary 
by the application of a positive velocity (i.e. from left to right) to the fluid. 
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Figure 5.2.: Change of wave-form calculated for (h — do M? p? > 0) 


we get 


dt — l 3. (40 M? p? — h) 


- (40 M? p? —h) -snpx-cnpx-dnpz. (55) 


. — 2. — 2 
oh = ph 7 [a — op? - (2 To pot 


Supposing then 
2: [a — op? - (2-—4M”’)| =3* (40M? *—h) ; 


we have 

dy 3qoph_ 

dt ss 
Here Figure 5.2 shows the change of form calculated for the case h — 4a M? p? > 0. 
When h — 4c M? p? = 0, the waves are stationary in accordance with Equation 3.7, 
whilst for h — 40 M? p? < 0 they become steeper behind; and this last result, since 
p is inversely proportional to the wave-length, may be stated by saying that cnoidal 
waves become less steep in front and steeper behind when, for a given modulus 
and a given height, their length is smaller than the one required for the stationary 
wave of this modulus and height. In proportion as M is taken smaller the cnoidal 
waves more and more resemble sinusoidal waves. They would take the sinusoidal 
form for M = 0, but then an infinitely small wavelength would be required for 
the stationary case. For this reason sinusoidal waves may always be considered as 
cnoidal waves whose length is too large to be stationary, that is, they are always 
becoming steeper in front. 


(40 M? p* — h) - sn® px - cn px - dn px (5.6) 
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Figure 5.3.: Change of wave-form calculated for sinusoidal waves 


5.3. Sinusoidal Waves 


This last result is easily verified by direct application of Equation 2.12 to the 
equation of a train of sinusoidal waves: 


20x 
=A-si —]|; 
" sin (7) , 


for, supposing 


ona 
a= yo 
we obtain 
dn — 3qo7A? ee Ame \ 
dt —- 21 N , 


and from this the change of form indicated in Figure 5.3 is easily calculated. 


5.4. More complicated Cases 


For the sake of curiosity, we represent by means of the following figures the change 
of form for some more complicated cases. 
Figure 5.4 and Figure 5.5 refer to the equation 


2 1 4 
n = A,-sin (=) + gt - sin (=) : 
Al 


In Figure 5.4, ++ is supposed to be small compared with (2), as is the case with 


waves of extremely small height. In Figure 5.5 we suppose Eig to be small in 
regard to a Generally for more complicated forms of waves these two cases have 
to be discriminated. When there is a moderate proportionality between the two 
fractions, the result is still more complicated. 
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Figure 5.4.: Sinusoidal wave — more complicated case 
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Figure 5.6.: Sinusoidal wave — more complicated case 
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5. Deformation of Non-Stationary waves 


Finally, Figure 5.6 refers to the equation 


2 1 4 
7 = A, - sin (=) — pel - sin (=) ; 
L 


in case that ( a is the smaller fraction. 
It is worthy of remark that all these waves grow steeper in front. 
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6. Calculation of the Fluid Motion 
for Stationary Waves to the 
Higher Order of Approximation 


In order to remove every doubt as to the existence of absolutely stationary waves, 
we will show how by development in rapidly convergent series the state of motion 
of the fluid belonging to such a wave-motion may be calculated. Expressing again 
the horizontal and vertical velocity of a particle by means of the series given in 
Equation 2.1 and Equation 2.2, which fulfil all the conditions for the interior of the 
fluid, we have only, neglecting capillarity, to satisfy the surface-conditions 


on 
Ox 


VUzy= U1 


(6.1) 
and 
us + v; + 2gn = const. (6.2) 


For the case of cnoidal waves, which is the general one, we have found as a first 
approximation 


(#) - an (h n)- (K+). 


But now, to obtain higher approximations, we assume, indicating by accents 
differentiation with respect to x: 


n? =an-(h—n)-(k+n):(1+bn+en*+---), (6.3) 


and 


f=qtrnt+sy +ty? +unt+--- (6.4) 
On writing out Equation 6.3, neglecting such terms as are of a higher order than 


the fourth compared with 7, h, and k, which latter quantities are of the same order, 
we obtain 


n” = ahkn + {a-(h—k) + abhk} -1? +{-a+ab-(h—k)}-7%—abn* (6.5) 


6. Calculation of the Fluid Motion ...to the Higher Order of Approximation 


and by differentiation 


1 
7! = salik + {a (hk) + abhk} + + {50+ sab. (h— m} 1° — 2aby’ . (6.6) 


From Equation 6.4, by successive differentiations and substitutions, retaining all 
terms up to the third and the 35% order, we deduce: 


f' = (r + 28n + 3én’) - 1; 


1 
"= sarhk + {ar (hk) + abrhk + 3ashk} «1 


3 3 
+f Sart Sabr (i =k) +4as-(h—K)} a+ {abr — Sas} 1 
f" = [ar -(h—k) + abrhk + 3ashk 
+ {-—8ar + 3abr - (h —k) + 8as- (h —k)}-1 + {—6abr — 15as} - 1 | ns 
; 1 
f? = 5a°rh- (h—k) + {er “(hk)” — sernich 1) 
— 2 ee -{h—k}-? + ye 
2 2 
v 2 2 39 2 49 9 9] i 
f? = |a°r-(h—k) — se rhk — lba°r- th —k} n+ Sarry 1; 


where 77’ is a quantity of the order of 3. 


Substituting these values in Equation 2.11, where y = 1+ 7, we have, retaining 
terms of the third order: 


_ 1 2 W " 1 4 au 
u=f 5 fi -ln-f t+ a4" f 
=q- part hk + ertihk -(h—k) 


+ " — care -(h—-—k)- Lae - Base hk — atht 
2 2 2 2 


Fe ae ae a ee 
52 (h—k) Tk hk? -n 


+ ‘3 + “aly — Sabri? -(h =k) —2asl’?-(h —k) 


—arl-(h—k) — 2 art -(h—- m} 9 


5 3 5 
| | 2 2 | 2,74 | 23 
! + abrl* 4 nasl sar t ig@ rl \ ie (6.7) 
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6. Calculation of the Fluid Motion ...to the Higher Order of Approximation 


We find in the same way, including terms of the 35m order: 


= ! / ls m 149 m tis v 
Se Sr Pr rl ae 


6 


29) 1 3. (h— pb) a abenp + dene — 2 ¢2rld- (hb — by? 
=H [rt + gant (h — k) abrl hk + sasl hk 720° 7 (h — k) 


Dope By ee ey ee ee 
+e50 8! hh} 2sl—r sari sobrl (h—k) 


+ Fast -(h—-—k)+ sar? -(h-—k)+ srl? -(h—- m} “1 


5 3 3 
{3H = 2s — abril? = 5st” _ sul — sere} : ta 


If now we write, in accordance with Equation 6.1 


m= SAL Bat Oy + Dit, 


we have from Equation 6.7 and Equation 6.8: 
A=q- 2 acl hike + + Prlthk -(h—k) 
4 48 
1 
=—rl+ gue -(h—k) 
a 1 br? hk + Past ie - See -(h—k)?+ > 25h k 
6 2 120 80 
a Deep is eos. 8 
BSF gar (h—k) sobrl hk Zasl hk sarlhk 
Lo dh B\2 2 2nd 
+ 54° rl*-(h—k) 16% eh 
1 1 
= —2sl—r — art + subrl? -(h—k) 
eee | 2 1 os 
C=s+ “alr - sabre -(h —k) — 2asl? - (h — k) 
—arl-(h—k)— 2 art -(h—k) 
= —3tl — 2s — abri® — Dep — ae — 3 o2yf5 
2 2 16 
5 


5 3 
D=t-+abri’ 4 5st” | sand Te 


(6.8) 


(6.9) 


(6.10) 


(6.11) 


(6.12) 


(6.13) 
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6. Calculation of the Fluid Motion ...to the Higher Order of Approximation 


Moreover, since Equation 6.2 may be written in the form 


ul: (1 + 77) +297 = (A + Bn + Cn? 4 Dn)* -(1+ ahkn 
t+a-(h—k)-7°- an’) + 2gn = const. 


we readily obtain 


2AB 4 
2AC + B? +a 


| ahkA? + 29g =0, 
-(h—-k)-A?=0, 
+ 2BC —aA? =0. 


2AD 4 


(6.14) 


(6.15) 
(6.16) 
(6.17) 


From Equation 6.10, Equation 6.11, Equation 6.12, Equation 6.15, Equation 6.16, 
Equation 6.17, the six quantities q, r, s, t, a, and b may be calculated, and if 
we had retained everywhere terms of one higher order, we might have got eight 
equations with eight unknown quantities, and so forth. 
By a first approximation we readily obtain from Equation 6.10 till Equation 6.13: 


psa 

l 

q 1 
$= 7 t gad 

q 1 
t= — 

B 3 
A=q 
hoes 

l 

q 1 
ye 2 oe 
_ a ,2 
=p ogee 


2 1 
— good ga al 


and then from Equation 6.15 till Equation 6.17: 
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6. Calculation of the Fluid Motion ...to the Higher Order of Approximation 


Proceeding to the second approximation, we find 


and then again from Equation 6.15 and Equation 6.16: 


eng (iu 050) 


Se) 


q (: | (h—k) 9 (h—-kP 98 *) 


an a 21 0 P80 - 
3 hk 

SE ges at 

pu_%42 (h—k) 21g (h—k)? 124 hk 

> sy oy l 20 1 [2 51 P 

>» (,,(h-k) 1 (h—k)? 38 hk 

v=9 (14 i 20° P 20°? 


By means of these results we may now readily obtain from Equation 2.1 and 
Equation 2.2 expressions for u and v including respectively the terms of the 2™4 


and 2 sn order. 


They are: 
(h-k) 3 (h—k)? 33 hk 
— . 14 a . — . 
v= val ( a 20 - 40° 2 
1. a))\ 2 Tet 
eg 1° 4P 
3 hk 3 (h—k)n 9 1 
| eta 7 1712) a) 
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6. Calculation of the Fluid Motion ...to the Higher Order of Approximation 


v=n #4(1 | ee TD yy (-“* | st) ath, (6.19) 


ae (1-3-9) een erm (143-2). 20 
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¢. Calculation of the Equation of the 
Surface 


We will now show how for the equation of the surface of a stationary train of waves 
a more correct expression than Equation 3.7 can be deduced. For this purpose, 
we have to integrate the differential equation Equation 6.3, or rather we have to 
prove that a series can be given which solves this equation to any desired degree of 
accuracy. Now such a series may be obtained in the following manner. Let 


1 h 
mM = h,-cn’ ae a: (hy + ky) [sean = a (7.1) 


represent the solution of an equation 


ny = am: (hi —m)- (ki +m) (72) 


where h, and k, have values which are slightly different from those of h and k in 
Equation 6.3; then these values and the coefficients a, 3, &c., of a series 


n=am + Bni+ynitéontt+-- (7.3) 


may be determined in such a way! that this series Equation 7.3 satisfies the equation 
Equation 6.3. 

Indeed, substituting Equation 7.3 in Equation 6.3 and taking into account Equa- 
tion 7.2, Equation 6.3 reduces to 


(a + 26m + 39n{ 4 +)? (Ay — m) + (ka +m) 


=(a+ 8m +m +-++)-(h-om — Bui — yn +++) 
-(k+am + Bnitynit---)- (1+ bam + (68 +ca) ni +--+), 


and it is only necessary to equalize the coefficients of the corresponding terms of 
both members of this equation. 


'The coefficient a in Equation 7.2 might also have been chosen slightly different in value from a 
in Equation 6.3, but this would only have introduced an unnecessary indeterminateness in the 
solution. 


7. Calculation of the Equation of the Surface 


If we retain all terms to the fourth order, we find in this way, after some reductions: 


a? - (hy — ky) —a?-(h—k)— (ba? — 38) e005) 
—a* + a* — (ba* — 2078) - (h — k) — (ca® — 2ba*B + 887 —5ay)-hk =0 (7.6) 
—4aB + 3078 + ba* — (ca? + 3ba*B — 36° — day) -(h-k)=0 (7.7) 
—48? — bay + ca* + 4ba® 8 + 308? + 3077 = 0. (7.8) 
To a first approximation these equations are satisfied by taking 

hi =h 

ky — k 
aad (7.9) 

B=b 

1 
Y= b? + ar ‘ 


If then we substitute in Equation 7.4, Equation 7.5, Equation 7.6, and Equation 7.7 


hy =hre 

kj, =kt+e 
a=lt+a, 
B=b+f,, 


where a and £ are quantities of the first, « and e of the second order, we find from 
these equations by second approximation: 


€ = —bhk 
ze = bhk 
Qgn= —b- (h a k) (7.10) 


1 
A= (-2" + 5°) -(h—k). 
Substituting as a third approximation: 


hy =h—-bhk+e 
ky =k+bhk+ x, 
a=1-b-(h-—k)+a2 
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7. Calculation of the Equation of the Surface 


we obtain finally: 


1 
c= gchk -(—h + 2k) 


1 
a2, = gchk - (2h — k) (7.11) 


2 
a2 = (08 — Se) (t= ne +2). 


Hence, the equation of the surface of the waves is, including all terms of the third 
order: 


n= fib. (h— b+ (Ze) «(82 — nk+ #2] mm 


1 1 
+ + (~20" + ) (n=) “m+ i + Fe] nee (7.12) 
where 
=) He a a: (hy +k) | modM = hh (7.1) 
m= 2 . : i hy t+ ky 
and 


1 
hy =h—bhk + =chk - (—h + 2k) 
: (7.13) 
ky = k-+ bak + schk- (2h— k). 


Here, according to Equation 6.20, 


(7.14) 
b=—+---, 


whereas the value of c and more correct expressions for a and 6 could only have been 
obtained by means of still more tedious calculations, which we have not executed. 
If we confine ourselves to that degree of approximation for which all the calculations 
have been effected, we may write for the equation of the wave-surface: 


y= [pS Ons oR (7.15) 


maheonts. (1- = =4)) .s. 3: (h+k) (7.16) 


ol 


7. Calculation of the Equation of the Surface 


m= (1-4). (7.17) 


For the solitary wave, when k = 0, we have? 


3h oe 
n= [1-3] m+ gen (7.18) 
7 ey ee ah 
m = h-sech 5 1 a}? iB (7.19) 


January 1895 


Cosme 5 


? Another close approximation of the surface-equation of this wave has been deduced by McCowan, 
[7], page 48. 
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A. The Korteweg—De Vries Equation 
in Modern Notation 


In modern notation the KORTEWEG-DE VRIES equation reads 
Ou+6ud,u+ O2u=0. (A.1) 


The KdV equation is a nonlinear, dispersive partial differential equation for a 
function wu of two dimensionless real variables, x and t, which are proportional to 
space and time respectively. 0, and 0; denote partial derivatives with respect to x 
and t. 

It can be derived from Equation 2.12 


OR 8. GO EL eg Te ORF 
at a fee (51 Rg ge ye et?) 


through the following substitutions: 


v= sé -t (A.2a) 
a (A.2b) 


Q. (A.2c) 


Differentiation of Equation A.2a with respect to Ot and Equation A.2b with respect 
to Ox yields 


ee ad (A.3a) 
OL = sue - Ox". (A.3b) 

Equation A.2c is rearranged to give 
n= 2u— a (A.4) 
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A. The Korteweg—De Vries Equation in Modern Notation 


From Equation A.3b it follows that 
Or = e+ On". (A.5) 


Inserting all this in Equation 2.12 yields 


a b= aVt ecm (FF) 


ae (2 - 5a) bees 2 eo) (A.6) 


3 3 a: Ox! 


Cancelling terms gives 


0 (2u — 2a) a (3 Ne 
ov ~C “ND (20 7 5°) 
Q 0? (2u — 2a) 
+2a- (20 — 5°) aD 3 (A.7) 


Differentiation of the left term in Equation A.7 gives: 


O (Qu — 2a) _ O(2u) 
Ot! Ot! 


O (—2a) ae OU 
ail 2 ai (A.8) 


Differentiation of the right term in brackets on the r.h.s. of Equation A.7 gives: 
0? (2u— Fa) — 0?(2u) | 0 (—$a) | ‘5 Oru 


12 = 2 «| 12 = 12° 
Ox Ox Ox Ox 


(A.9) 


The two remaining terms in brackets on the r.h.s. of Equation A.7 are multiplied 
to give 


3 2\7 2 3 8 4 
a (20 = 5°) 20h (20 - 5°) 5 G — gua oe 5°”) 


4 
4 Los es 2 
+ 4ua 3 
2 4 
= 6u? — dua + —a? + 4ua — —a? 
3 3 
2 
= 6u? — oe (A.10) 
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A. The Korteweg—De Vries Equation in Modern Notation 


Inserting Equation A.8, Equation A.9 and Equation A.10 in Equation A.7 yields 


du OA ae Oru 
2. a= pe (6 =30 +255] 


and differentiating again gives 


Ou 0 (6u7) | 0 (—2a7) Pu 


gpg ge Es 
Ou Ou OPu 
rn a oe 


Multiplying both sides of Equation A.12 with 5 finally gives 


Ot! Ox! = Ox 
as a du 8 
U U U 

ap be Bye 


(A.11) 


(A.12) 


(A.13) 


(A.1) 
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B. Jacobi Elliptic Functions and 
Elliptic Integrals 


The JACOBI elliptic functions! are standard forms of elliptic functions. The 
three basic functions are denoted cn(u, k) (cosinus amplitudinis), dn(u, k) (delta 
amplitudinis), and sn(u,k) (sinus amplitudinis)*, where k is known as the elliptic 
modulus. From these three basic functions the remaining 9 functions are derived. 
In sum there are 12 JACOBI elliptic functions. 

The name elliptic function or integral stems from the fact that the perimeter of 
an ellipse is most simply expressed in terms of the complete LEGENDRE elliptic 
integral of the second kind, and this is why the functions are so called. 


J R(w, x) dz and f R(./ f(x), x) dz respectively are a general form of an elliptic 
integral, where R(w, 2x) is a rational function of x and w, w? is a function of z, 
ie. f(x) = w? or / f(x) = w respectively and f(x) is a polynomial function 
that is cubic or quartic in x, and R(w, x) contains at least one odd power of w, 
and w? has no repeated factors. Furthermore elliptic integrals can be viewed 
as generalizations of the inverse trigonometric functions, known as elliptic 
functions, but these cannot, in general, be expressed in terms of elementary 
functions. 

However, with the appropriate reduction formula, every elliptic integral 
can be brought into a form that involves integrals over rational functions 
P(w, x), Q(w,x) and the three LEGENDRE canonical forms, i.e. the incomplete 
elliptic integrals of the first, second and third kind denoted by F'(¢,k), E'(¢,k), 
and II(n; ¢,k), respectively. In case of ¢ = 7/2, then the integrals are called 
complete elliptic integrals and are denoted K(k), E(k), and Il(n;k). 


JACOBI elliptic functions arise from the inversion of the incomplete elliptic 
integral u of the first kind (type I), 


(B.1) 


uU= 


F( y= f° = 
0 /1 —k?- sin2t 


‘Taken mainly from https://mathworld.wolfram.com/, https://dlmf.nist.gov/22 and Wikipedia 

?This notation was introduced by CHRISTOPH GUDERMANN (* 25. March 1798, Vienenburg; 
{ 25. September 1851, Miinster) in 1838. The terminology stems from the old, and outdated, 
notation of JACOBI: snu = sinamu, cnu = cosamu, dnu = Aamu 


B. Jacobi Elliptic Functions and Elliptic Integrals 


where 0 < k? <1, k = mod u is the elliptic modulus, and ¢ = am(u, k) = am(uw), 
with —1/2 < @ < 7/2 is the JACOBI amplitude, an infinitely-valued function of u, 
giving 
1 
o= plus Ry = amd &) (B.2) 
From Equation B.2 the relationships 


sin(d) = sin(am(u, k)) 


= snd; 4) 
cos(@) = cos(am(u, k)) 
= cn(u, k) 
V1 — k?-sin?(¢) = V1 —k?- sin?(am(u, k)) 
= dn (uk). 


follow. 
Equation B.1 is the trigonometric form of the elliptic integral type I; substituting 
6 =sint and x = sing, one obtains the LEGENDRE normal form: 


id dé 
oy, Ja-®) 0 — Pe) 


The complete elliptic integral of the first kind is given by K(k) = F (5, k), or 


-{ dt -[ da 
0 V1—k?-sin?t Jo \/(1— 6?) - (1 — 6?) ’ 


if one substitutes x = sin(7/2) = 1 and sint = 0. 

In general the integrals with a lower integral limit of 0 are called incomplete elliptic 
integrals. Is also the upper integral limit >, in the case of types I and II we speak 
of complete elliptic integrals. 


Incomplete elliptic integrals of the second kind F (type IJ) in trigonometric form 


are defined as 
b= i V1—k?-sin*tdt. (B.5) 
0 


Substituting 0 = sint and x = sin @, one obtains the LEGENDRE normal form: 


Els, k) Nee vi-k@ 


a (B.6) 


3sint = 0, > do cost, and therefore dt = 


dt 
dé 
V1—62° 


cos a which gives, recognizing sin” t + cos?t = 1 
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B. Jacobi Elliptic Functions and Elliptic Integrals 


Complete elliptic integrals of the second kind F in trigonometric and LEGENDRE 
normal form are defined as 


y= [ ake rdee [oS asl (B.7) 


(/1— 


with sin(7/2) = 1 and 6 = sint. 


Incomplete elliptic integrals of the third kind II in trigonometric form are defined 


co) 
Ce oe / : a =) (B.8) 


1 —nsin’ t) /1-—(sintsina)? 


as 


The number n is called the characteristic and can take on any value, independently 
of the other arguments. Substituting 0 = sint and x = sin@, one obtains the 
LEGENDRE normal form: 


Se a dd 
Geis i (1 — nO?) \/(1 — k262) - (1 — 62) ey 


Complete elliptic integrals of the third kind II in trigonometric and LEGENDRE 
normal form are defined as 


$/2 1 dt 

as - © ~ / (1 —nsin? t) \/1 — k2- sin? t em) 
tt dé 

OA i (1 —n6?) \/( — k62) - (1 — 6) Bat) 


with sin(7/2) = 1 and 6 = sint. 
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B. Jacobi Elliptic Functions and Elliptic Integrals 


The functions sn(u, &), cn(u,&), and dn(u, &) are doubly periodic generalizations 
of the trigonometric functions satisfying, for the limiting cases k = 0 and k = 1, 
the following relations: 


function k=0 k=1 


sn(u,k) sin(u) tanh(u) 


cn(u, k) cos(u) sech(w) = ae 


dn(u,k) 1 sech(u) = can 


Table B.1.: Limiting forms of JACOBIAN elliptic functions with regard to the elliptic 
modulus k 


I.e., When k = 0 or k = 1 the JACOBI elliptic functions degenerate into trigono- 
metric or hyperbolic functions. 


The three main JACOBI elliptic functions sn(u), cn(u), and dn(u), from which 
all others are derived are periodic in K(k) and K'(k). 
In general, the double periodicity of the JACOBI elliptic functions may be expressed 
as: 

pq(u + 20K (k*) + 2i8K (1 — k?) , k’) = (-1)’paq(u, k”) (B.12) 

where a and ( are any pair of integers. K(-) is the complete elliptic integral of the 
first kind, also known as the quarter period. The power of negative unity (7) is 
given in the following table: 


q 
c s n 
c 0 B a+ B a 
s B 0 a a+ B 
. nla+ 6 ay 0 B 
d a a+ 6 B 0 


Table B.2.: Values for the power 7 of negative unity 


When the factor (—1)7 is equal to —1, the equation expresses quasi-periodicity. 
When it is equal to unity, it expresses full periodicity. It can be seen, for example, 
that for the entries containing only a when a is even, full periodicity is expressed by 
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B. Jacobi Elliptic Functions and Elliptic Integrals 


the above equation, and the function has full periods of 4K (k?) and 2i K(1 — k?). 
Likewise, functions with entries containing only 6 have full periods of 2K (k?) and 
4i K(1 — k?), while those with a + 6 have full periods of 4K (k*) and 41 K(1—k?). 
For the three main JACOBI elliptic functions, along with their zeros and poles, the 
result is shown in Table B.3. 


function periods ZeTO pole 

sn(u,k) 4h,2i Kk" 2B In 2mK + (2n+1)ik’ 

en(u,k) 4K,2(K+iK') (Qm4+1)K4+2nik’ 2mK + (2n+ 1)ik’ 

dn(u,k) 2h,4i Kk" (Qm4+1)K+(Qn4+1)iK’ ImK+ (2n+1)ik’ 
n,mEeZ 


Table B.3.: Periodicity, zeros and poles of sn(u), en(w), and dn(w) 


K(k) is the complete elliptic integral of the first kind as stated above with 
K'(k) = K(k’), k' = V1 —k? (the complementary modulus) and therefore K(k) 
and K'(k) are interrelated through the elliptic integrals 


n/2 n/2 
222 f ie i: aS 


0 1—k?2 sin? t 


dt 
1—k?)-sin?t 


E.g. sn(u,k) has zeros at u = 0 and u = 2K and the poles at u = ik’ and 
u=2K +iK’. 


The cnz, dna, and sn functions may also be defined as solutions to the 
second-order nonlinear ordinary differential equations 


d? 

a = 2B) yw 

2 

Co = — (1-28) y — 2h?98 
dy 2 2,3 


respectively. 
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B. Jacobi Elliptic Functions and Elliptic Integrals 


The standard JACOBI elliptic functions satisfy the identities 
sn?u+cn?u=1 
k*sn?u+dn?u=1 
kon? ut k? = dn? u 
cw u+k?sn?u =dn*u. 


Special values include 


en(O.k) =cn(0) 1 

enA(A).b) =n Kk (b))=0 
dn(0; &):=dn(0): = 1 

dn(K(k),k) = dn(K(k)) = k! = V1 — k? 
sn(0, 4) =sn(0) =1 

sn(K(h).k) = snk (kh) ):=0 


where K = K(k) is a complete elliptic integral of the first kind and k’ = V1 — k? 
is the complementary elliptic modulus, and 


cn(u, 1) = sech u (B.13) 
dn(u, 1) = sech u (B.14) 
sn(u, 1) = tanhu, (B.15) 


as stated above. 


Al 


B. Jacobi Elliptic Functions and Elliptic Integrals 


In terms of integrals, 


< 

I 
a 
eB 
s 
a 


Gay Garr a 
(Ga @CS)) oat 


(1 — #2)" (h? + 2?) at 


I| 
a ee 


(Psa) OP) ae 


a 


(1 a 7); (t? = i ea dt 


(Paty (eRe ys ai 


D 
ie) 
_~ 
s 
= 
YS 


(iar) eee) a 


I 
pSyos eS eS 


fe) 
n 
—~ 
s 
oy 
ale 


(t? + as (i? ifs ea ee dt 


sd(u,k) 


(1 _ Rae (1 4 ey? dt 


Cay eae a 


cd(u,k) 
Gee) Osh) a 


BR 


By go 


(a2 ak) ds 
de(u,k) 


JACOBI elliptic functions addition formulas include (where, for example, sn(u, k) is 
written as sn u for conciseness) 


snucnudnvy+snvucnudnu 


Saat) 1—k*sn2 usn2v 
_ chucny — snusnudnudnv 
Se ie 1—k?sn2usn?2 v 
dnudnv — k? snusnvcnucnv 
dn(u+v) = 


1—k?sn2 usn2v 
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B. Jacobi Elliptic Functions and Elliptic Integrals 


Extended to integral periods, 


sn(u+ K) = arr 
cn(u+ Kk) =— cent 
k! 
dn(u+ Kk) = Aer 
sn(u+ 2K) =—snu 
cn(u+ 2k) =—cnu 
dn(u+ 2k) =dnu 


For complex arguments, 
Paar sn(u, k&) dn(v, k’) 
~ 1 —dn?(u, k) sn?2(v, k’) 
. (a k) dn(u, k) sn(v, k’) en(v, —)) 
| 1 — dn?(u, k) sn?(v, k’) 
.. en(u, k) en(v, k’) 
OS 1 — dn?(u, k) sn?(v, k’) 
att (= k) dn(u, k) sn(v, k’) dn(v, “)) 
1 — dn?(u, k) sn?(v, k’) 
/ / 2 / 
re ee dn(u, Dea k') dn(v, k’) (* sn(u, f) en(u, k) sn(v, k ‘) 
1 — dn*(u, k) sn?(v, k’) 1 — dn*(u, k) sn?(v, k’) 


Derivatives of the JACOBI elliptic functions include 


d(sn.u) =cnudnu 
du 

ox) = —snudnu 
du 

adam) = —k’snucnu 
du 


Double-period formulas involving the JACOBI elliptic functions include 


2snucnudnu 


sn(2u) = 


1—k?sntu 
(2u) 1—2sn?u+ k?sn*u 
n — 
peer 1—k?sntu 
1 — 2k? snu + k? sn*u 
dn(2u) = 
n(2u) 1—k’%sntu 
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B. Jacobi Elliptic Functions and Elliptic Integrals 


Half-period formulas involving the JACOBI elliptic functions include 


dn (5K ) — Vk 
Squared formulas include 
aren 
cats OT aa 
an ay 


TAYLOR series of the JACOBI elliptic functions were considered by HERMITE (1863), 
SCHETT, [11], [10] and Dumonr [5], 


1 1 

on(w, k) =1— Su" + 57 - (1+ 4k") ul — Se - (1 + 44k? + 16h!) uh + --- 
1 1 4y,4_ 1 2 4, 16\,6 
1 1 

sn(u,k) =u——- (1+?) wu mo + 14k? + k*) uh +... 
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